Anisotropy driven ultrafast nanocluster burrowing 
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We explore the occurrence of low-energy and low-temperature transient cluster burrowing lead- 
ing to intact cluster inclusions. In particular, the anomalously fast (ballistic) Pt nanocluster im- 
plantation into Al and Ti substrates has been found by molecular dynamics simulations using a 
tight-binding many-body potential with the 1 — 5 eV/atom low impact energy. Similar behavior 
has also been found for many other cluster/substrate couples such as Cu/Al and Ni/Ti, Co/Ti, etc. 
In particular, in Ni/Ti at already ~ 0.5 eV/atom impact energy burrowing takes place. At this 
few eV/atom low impact energy regime instead of the expected stopping at the surface we find the 
propagation of the cluster through a thin Al slab as thick as ~ 50 A with a nearly constant speed 
(oc 1 eV/atom). Hence the cluster moves far beyond the range of the impact energy which suggests 
that the mechanism of cluster burrowing can not be explained simply by collisional cascade effects. 
In the couples with reversed succession (e.g. Ti/Pt, Al/Pt) no burrowing has been found, the clus- 
ters remain on the surface (the asymmetry of burrowing). We argue that cluster penetration occurs 
at few eV/atom impact energy when the cluster/substrate interaction is size-mismatched and mass 
anisotropic atomically. 
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I. INTRODUCTION 

Cluster deposition on solid surfaces has been the 
subject of intense atomistic level studies in the last 
decades 1-7 . The reason of the considerable interest is 
due to that recent progress in this field offers a great 
technological potential for the application of clusters in 
various areas, such as the smoothening and cleaning of 
surfaces 8 , the improvement of magnetic properties 9 , size- 
dependent catalytic activity 10 , thin film growth with a 
well-controlled grain size 1-3 or the formation of cluster 
assembled nanowires 11 . 

Cluster-surface interaction as well as the mobil- 
ity of clusters on surfaces have also been intensively 
studied 3 ' 12-15 . Transient mobility (TM) of clusters has 
been known for a while on surfaces 2 ' 12,17 ' 18 but not in 
the bulk. Hereby we report on the possible occurrence 
of TM of various clusters perpendicular to the surface 
leading to transient cluster burrowing into various sub- 
strates. Until now only slow cluster burrowing with soft 
landing has been studied in few systems which takes 
place at least in a nanosecond time scale well above room 
temperature 13-15 . The high energy ion-bombardment in- 
duced burrowing of Pt nanoparticles into Si02 has also 
been reported recently 16 . 

Cluster impact phenomena have been intensively stud- 
ied by computer atomistic simulations to understand the 
processes which govern the creation of nanoscale struc- 
tures on the surface 1-4,6 ' 7,13-15 . However, at best of our 
knowledge no low-energy (few eV/atom or less) and low 
temperature cluster beam deposition with transient im- 
plantation rate has been reported until now in which the 
cluster remains intact in the bulk. Implantation of clus- 
ters with few tens of A depth usually requires at least 



few tens of eV impact energy 2 or even higher energy is 
required for deeper penetration depths 6,7 . Cluster im- 
plantation at much lower energies could offer great tech- 
nological possibilites for the production of nanostructures 
due to the much less destructive conditions. 



According to our recent findings that single atomic de- 
position in Pt/Al(lll) 19 and interfacial intermixing in 
Pt/Ti bilayer 20 are anomalous in many respect, e.g. tran- 
sient inter-layer atomic transport rates have been found, 
it could be that cluster transport in systems with the 
same atomic constituents could also be anomalously fast. 
Moreover, it has also been shown recently that in mass- 
anisotropic bilayers the enhancement of interfacial broad- 
ening and intermixing takes place 22 ' 23 and it could also 
be that mass-anisotropy does influence the propagation 
of deposited clusters. 



In this article we report on the occurrence of ultrafast 
cluster transport (burrowing) in Al (and in Ti, Ag and 
in Au) bulk upon low-energy impact (< few eV/atom) at 
~ K predicted by atomistic molecular dynamics (MD) 
simulations. The employed interpolated crosspotentials 
have been fitted to ab initio calculations. We show that 
the Pt cluster initialized with few eV/atom impact en- 
ergy can move ballistically in the bulk of Al until few 
tens of ps. In other deposition events with various clus- 
ter/substrate couples we find similar behavior when the 
cluster to substrate interaction is atomically size mis- 
matched and/or mass anisotropic in such a way that the 
cluster is composed of heavier and smaller atoms. 
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II. THE SIMULATION METHOD 
A. General properties 

Classical constant volume tight-binding molecular dy- 
namics simulations 26 were used to simulate soft landing 
and low-energy cluster impact of various nanoclusters on 
(111) surfaces of Al, Ti, Ag and Au at ~ K initial 
temperature using the PARC AS code 4 ' 13 - 27 . The MD 
code has widely been used for the study of various atomic 
transport phenomena in the last few years 4,13 ' 19,20 ' 22,27 . 
Although we carry out simulations at ~ K, we find 
a substantial local heating up in a local surface region 
of Al, hence the correct dissipation of the emerged heat 
should be handled using temperature control. A vari- 
able timestep and the Bcrcndsen temperature control 
is used at the cell border 28 . The simulation uses the 
Gear's predictor-corrector algorithm to calculate atomic 
trajectories 27 . The maximum time step of 0.05 fs is 
used during the operation of the multiple time step algo- 
rithm. The system couples to a heat bath via the damp- 
ing constant to maintain constant temperature condi- 
tions and the thermal equilibrium of the entire system 28 . 
The time constant for temperature control is chosen to 
be t = 70 fs, where r is a characteristic relaxation 
time to be adjusted 28 . The Bcrcndsen temperature con- 
trol has successfully been used for nonequilibrium sys- 
tems, such as occur during ion-bombardment of various 
materials 4 ' 13 ' 20 - 22 ' 27 . 

For simulating deposition it is appropriate to use tem- 
perature control at the cell borders. This is because it is 
physically correct that potential energy becomes kinetic 
energy on impact, i.e. heats the lattice. This heating 
should be allowed to dissipate naturally, which means 
temperature control should not be used at the impact 
point. Periodic boundary conditions are imposed later- 
arily. The observed anomalous transport processes are 
also observed without periodic boundary conditions and 
Berendsen temperature control. Further details are given 
in 27 and details specific to the current and similar sys- 
tems in recent communications 19,20 ' 22 . 

The top of the simulation cell is left free (the free sur- 
face) for the deposition of Pt atoms. The bottom layers 
are held fixed in order to avoid the rotation of the cell. 
Since the z direction is open, rotation could start around 
the z axis. The bottom layer fixation is also required to 
prevent the translation of the cell. In few cases, however, 
we do not fix bottom layers, when the penetration of the 
deposited cluster through thin substrate films has been 
demonstrated. No rotation of the cell has been found in 
these cases. 

The size of the typical simulation cell for the Pt / Al sys- 
tem is 80 x 80 x 42 A 3 including 16128 substrate atoms 
(with a fee lattice, 15 active MLs are supported on 3 
fixed bottom monolayers (MLs)). We also tried sub- 
strates with different lateral sizes (80 - 160 A ) and slab 
thickness (42 — 100 A ) up to ~ 150000 number of atoms. 



TABLE I: The parameters used in the interpolated tight bind- 
ing potential (TB-SMA) given in Eqs . (l)-(2) 26 The parame- 
ters of the crosspotential have been obtained as follows using 
an interpolation scheme 22 : For the preexponentials £ and A 
we used the harmonic mean A^ipt = (Aai x Apt) 1 ^ 2 for q 
and p we use the geometrical averages: qAiPt = (qai +qpt)/2. 
The first neighbor distance of the Al-Pt potential is given also 
as a geometrical mean of ro = (rjf* + r f f i )/2. In order to get 
better agreement between ab initio and semiempirical poten- 
tials parameters £ and ro have been optimized. In the case of 
Co/Au parameter q has also been fitted. 

£ q A 



P 



Al-Pt 2.7 3.258 0.191 9.612 3.0 
Ti-Pt 4.2 2.822 0.149 11.015 2.87 
Ni-Ti 0.8 1.416 0.052 14.209 2.72 
Cu-Al 2.1 2.397 0.102 9.785 2.85 
Co-Au 2.5 2.7 0.140 10.917 2.85 



Other cluster/substrates couples (Pt/Ti, Co/Au, Cu/Al, 
Ni/Ti) together with reversed succession (Ti/Pt, Ti/Ni, 
etc.) have also been constructed with similar system size. 
We tested cluster impact on much larger substrates and 
find no dependence of the anomalous atomic transport 
properties of the deposited atoms on the finite size of 
the simulation cell. The cubic cluster includes roughly 
<~ 1500 atoms (the length of cube edges: ~ 30 A). We 
also tested the penetration of various clusters in different 
substrates with similar lateral size and slab thickness and 
cluster size. The kinetic energy of the deposited particles 
is in the range of ~ 1 — 10 eV/atom (the velocity of the 
clusters normal to the surface iswsi 20-200A/ps). The 
cluster is initialized at <~ 8 A above the (111) surface of 
the substrate. 



B. The interatomic potentials 

We use the many-body Clcri-Rosato (CR) parametriz- 
ation of the tight-binding second-moment approximation 
(TB-SMA) interaction potential to describe interatomic 
interactions 26 . 

Within the TB-SMA, the band energy (the attractive 
part of the potential) reads, 



K = - 



j,rij<r c 



£ exp 



-2q 



ro 



- 1 



1/2 



(1) 



where r c is the cutoff radius of the interaction and ro is 
the first neighbor distance (atomic size parameter). 

The repulsive term is a Born-Mayer type phenomcno- 
logical core-repulsion term: 



El = A ^ ex P 



-P 



ro 



(2) 
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FIG. 1: The crosspotential energy (eV) for the dimers Al-Pt, 
Ti-Pt (inset Fig. la), Ni-Ti, Cu-Al (inset Fig. 2a) and Co/Au 
as a function of the interatomic distance (A) obtained by the 
ab initio PBE/DFT method. For comparison the fitted inter- 
polated semiempirical pote ntials (TB-SMA) are also shown. 
The fit procedure has been carried out by varying parameters 
£ and ro in order to get the best macth of the potential energy 
curves with the ab initio one. 



The parameters (£, q. A,p, r ) are fitted to experimental 
values of the cohesive energy, the lattice parameter, the 
bulk modulus and the elastic constants en, c\2 and C44 26 
and which are given in Table 1. The summation over j is 
extended up to fifth neighbors for fee structures 26 . The 
cutoff radius r c is taken as the third neighbor distance for 
all the interactions. We tested the Al-Al and the Al-Pt 
potential at cutoff radius with larger neighbor distances 
and found no considerable change in the results. This 
type of a potential gives a very good description of lattice 
vacancies, including migration properties and a reason- 
able description of solid surfaces and melting 26 . We also 
tested the Al-Al and the Al-Pt potential at cutoff radius 
with larger neighbor distances and found no considerable 
change in the results 23 . The cutoff radius is taken as the 
third-forth neighbor distance (r c w 10—15 A) for all the 
interactions which we find sufficiently large enough. 

For the crosspotential of substrate atoms and Pt we 
employ an interpolation scheme 22 using the geometrical 
mean of the elemental energy constants and the har- 
monic mean for the screening length are taken as in 
refs. 22 . The CR elemental potentials and the interpo- 
lation scheme for heteronuclear interactions have widely 
been used for MD simulations 14,21,22 . Recently a CR in- 
terpolated crosspotential has also successfully been used 
for Ti/Pt in agreement with our experimental results 20 . 
The scaling factor ro (the heteronuclear first neighbor- 
distance) is calculated as the average of the elemental 
first neighbor distances. The AlPt potential has been fit- 
ted to the measured effective heat of mixing in the cubic 



AlPt (AH w -100 kJ/mol) 22 . The melting temperature 
of 1870 K 30 is reproduced by our Cleri-Rosato crosspo- 
tential within the range of 1800 ± 100 K. In order to 
adjust AH in the Al-Pt potential (which is proportional 
to the strength of the interaction and to the heat of al- 
loying in the AlPt alloy) the preexponential parameter £ 
in Eq (3) is set to £ w 3.0 23 . Further details are given 
in 27 and details specific to the current system in recent 
communications 20,22,23 . 

In order to check the accuracy of the employed interpo- 
lated crosspotentials, the crosspotential energy has also 
been calculated for few of heteronuclear pairs (Al-Pt, Ti- 
Pt, Ni-Ti, Cu-Al, Co-Au) using ab initio local spin den- 
sity functional calculations 31 together with a quadratic 
convergence self-consistent field method. The G03 code 
is well suited for molecular calculations, hence it can be 
used for checking pair-potentials. The interatomic poten- 
tial V(dr) between two atoms is defined as the difference 
of total energy at an interatomic separation dr and the 
total energy of the isolated atoms 

V(dr) = E(dr) - E(oo). (3) 

The Kohn-Sham equations (based on DFT) 32 are solved 
in an atom centered Gaussian basis set and the core elec- 
trons are described by effective core potentials (using the 
LANL2DZ basis set) 33 and we used the Perwed-Burke- 
Ernzerhof (PBE) gradient corrected exchange-correlation 
potential 34 . First principles calculations based on density 
functional theory (DFT) have been applied in various 
fields in the last few years 35 . The obtained profiles are 
plotted in Fig. 1 together with our interpolated semiem- 
pirical many-body TB-SMA potentials for the various 
dimers. Our interpolated potential is partly fitted to the 
ab initio one: parameters £ and ro have been varied to get 
the best matching. We find that our fitted interpolated 
TB-SMA potentials match reasonably well the ab ini- 
tio one hence we are convinced that the TB-SMA model 
accurately describes the heteronuclear interaction in the 
Al-Pt dimer. Even if we use the interpolated TB-SMA 
without fitting, we get nearly the same results which in- 
dicates that the interpolation scheme is effective. We 
assume that the fitted diatomic potential is transferable 
for those cases when a cluster atom is embedded in the 
substrate. Heteronuclear interactions can be reasonably 
well described by dimer potentials during interfacial in- 
teractions while pair potentials are less accurate in al- 
loy phases where the local number of crossinteractions is 
much larger than 1. 

In Figs la-lc we show the ab initio DFT potential en- 
ergy curves and the fitted interpolated ones. In most of 
the cases we could reach a fairly good agreement which 
might increase the credibility of the employed approach. 
We also carried out simulations for couple of other clus- 
ter/substrate pairs and for which we use the simple inter- 
polation scheme for the crosspotential. Results for these 
pairs are shown in the discussion section. 
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FIG. 2: The cartoons of the simulation cells obtained after 
the impact of Co cluster on Au(lll) and Pt nanocluster on 
Al(lll). The atoms with lighter color are the cluster atoms. 
From right to left: fig. 2a: the top view of the simulation 
cell of Co cluster on Au(lll) impact event after 5 eV/atom 
impact, fig. 2b: the burrowed (embedded) Pt cluster in Al 
(crossectional view, leV/atom impact ), fig. 2c: the top view 
of the burrowed Pt cluster, fig. 2d: the crossectional slab (cut 
in the middle of the cell) of the Al(lll) substrate after the 
penetration of the implanted Pt cluster (5 eV/atom impact 
energy). Figs. 2a- 2c have been plotted after few tens of ps 
the cluster impact. In Fig. 2d the snapshot is plotted at 5 ps. 



III. RESULTS 

The low-energy deposition (1 — 5 eV/atom) of nan- 
oclusters of Pt on Al(lll) and Ti leads to the unexpected 
ultrafast cluster sinking and burrowing within few ps 29 . 
The process can be considered as an athermal transient 
transport which takes place in nearly K simulations 
within a ps at few eV/atom imapct energy. We find a 
similar behavior for other cluster/substrate couples such 
as Cu/Al, Ni/Al, Ni/Ti or Co/Ti, etc. Common fea- 
ture of these couples is the large atomic size or mass 
anisotropy (the host atoms are large and light and the 
cluster atoms are smaller and heavy). 

Normally, most of the clusters should stop and fall 
apart and spread on the surface or in the upper layers 
at this energy regime (see e.g. in the review article 2 or in 
refs. 3,6 ' 7 ). Actually, this feature of cluster impact can be 
used for thin film growth 1-3 . In this article, this has been 
demonstrated for the Co cluster deposition on Au(lll) 
in Fig. 2a and can also be seen at a web page 29 . In this 
case we find the fragmentation of the Co cluster upon im- 
pact on Au(lll) at 5 eV/atom impact energy. We expect 
for most of the low-energy cluster deposition events with 
various substrates and clusters the stopping of the clus- 
ter on the surface (with only weak penetration). Even 
for higher impact energies we find no penetration of Co 
to Au (up to few tens of eV/atom). At 5 eV/atom en- 
ergy the Co nanoparticle falls apart and an intermixed 



nanodot remains on the surface. 

Also, no burrowing occurs for the couples with reversed 
succession: e.g. the energetic deposition of Ti cluster on 
Pt substrate leads to trapping on the surface with some 
fragmentation and/or spreading on the surface (similar 
situation takes place for Al/Pt and Ti/Ni, etc.) 29 . Hence 
the ultrafast burrowing of clusters is asymmetric with 
respect to the interchange of the atomic constituents. 

Clearly we find a different and anomalous situation 
for Pt cluster on Al(lll) and for the other couples un- 
der study (Pt/Ti, Ni/Ti, Cu/Al). Instead of the normal 
cluster low-energy impact behavior the Pt cluster keeps 
its integrity during the impact and starts to burrow (pen- 
etrate) ballistically below the surface instead of stopping 
at the surface. Below 1 eV/atom the stopping or the par- 
tial penetration of the Pt clusters has been found which 
implies that the burrowing has a potential energy barrier 
of ~ 1 eV/atom. 

Penetration of clusters upon impact has only been 
found until now at much higher impact energy above 
~ 15 eV/atom which takes place together with the frag- 
mentation of the cluster 2 ' 6 - . Hence no burrowing pro- 
cess takes place and the higher implantation energy of 
the clusters lead to the serious damage of the cluster and 
the substrate. E.g. the amorphization of size-selected 
Au, Ag, and Si clusters have been found in graphite at 
few keV implantation energy 6,7 . Our finding is also dif- 
ferent from the cases found for Co cluster burrowing in 
Cu. In this case the cluster penetration (CP) is relatively 
slow process and takes place in a ns time scale and is a 
thermally activated process 13 . 

It should also be stressed that the Pt cluster remains 
also intact at higher impact energies. In Figs. 2b-2c the 
cartoons of the burrowing event is shown for the case 
when the cluster is initialized by 1 eV/atom kinetic en- 
ergy. The sinking of the cluster takes place in few ps 
(ballistic process). Increasing the impact energy up to 
5 eV/atom the sinking process continues and the cluster 
penetrates through a ~ 50 A thick Al or Ti slab (Fig. 
2d) in 10 ps. The cluster travels through the thin sub- 
strate slab with a nearly constant speed (~ 1 eV/atom) 
and remains nearly intact when leaves the bottom of the 
thin slab. The corresponding animations can be seen in 
a web page 29 . In Fig. 2d we show a snapshot of cluster 
implantation at 5 ps. 

Hence the stopping of the Pt cluster does not occur as 
it should be in most of other cluster impact events (such 
in the case of Co cluster on Au or for Al cluster on Al). 
In the case of Co cluster on Au(lll) we see the complete 
disintegration of the Co cluster on the surface (Fig. 2a) 
even at 5 eV/atom impact energy. At 5 eV/atom impact 
energy the Pt cluster penetrates far beyond the range of 
the deposited energy (that is few A) into the substrate by 
an order of magnitude larger implantation depth. There- 
fore, simple collisional cascade effects can not explain the 
enhancement of the vertical cluster mobility. Therefore, 
there must be a specific mechanism which weakens stop- 
ping effects for the Pt/Al, Pt/Ti for other couples. 
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FIG. 3: The average kinetic energy (eV/atom) of the pene- 
trating Pt cluster as a function of the simulation time (ps) 
during complete penetration through the Al slab (the slab 
thickness is ~ 50 A). Inset: The crossectional view of 
the crossectional slab of the burrowed Pt cluster in Al at 1 
eV/atom impact energy cut in the middle of the simulation 
cell. 



FIG. 4: The penetration depth (d in A) of various cluster 
impact events as a function of the impact kinetic energy. Er- 
ror bars denote standard deviations obtained during various 
events (statistics). Inset: The average kinetic energy of pen- 
etration of cluster atoms (eV/atom) through the substrate as 
a function of the impact kinetic energy (eV/atom). 



IV. DISCUSSION 



In Fig. 3 the average (downward) kinetic energy of 
the moving (burrowing) cluster per atom is also shown 
for a typical transient cluster burrowing (TCB) event at 
5 eV/atom initial kinetic energy. The cluster slows down 
to ~ 1 eV/atom after the impact and continues moving 
through the Al slab. Contrary to this, the Co cluster 
stops immediately after the impact and the kinetic energy 
of the cluster drops sharply to zero at the same initial 
impact energy. In the case of 1 eV/atom impact energy 
we find burrowing for Pt on Al and no further penetration 
is seen which indicates that complete cluster penetration 
has a barrier of few eV/atom. The crossectional view of 
the embedded cluster is shown in the inset Fig. 3 at 
the end of the simulation which shows us that the cluster 
remains largely intact during burrowing. Also, even no 
partial cluster burrowing occurs below ~ 0.5 eV/atom 
impact energy. 

In Fig 4 we show the calculated penetration depth d p 
of clusters obtained during various impact events mak- 
ing statistics as a function of the impact energy up to 10 
eV/atom. We find the nearly linear increase of d p with 
Eimp- This is in accordance with the finding of Praton- 
tep et al. 6,7 as they found the nearly linear scaling of 
the implantation depth of Au, Ag, and Si clusters into 
graphite as a function of the impact energy in the keV 
impact energy regime. We find this scaling relation at 
much lower impact energy (few eV/atom regime). In the 
Inset Fig 4 the linear scaling of the kinetic energy of 
penetration (E pen , propagation) can be seen as a func- 
tion Ei mp . These features show us that both d p and E pen 
scales linearly with Ei mp . 



Interestingly, the observed anomalous processes (tran- 
sient burrowing, the asymmetry of CP) are only weakly 
sensitive to the choice of the heteronuclear potential. We 
get similar results for the interpolated and for the fit- 
ted potentials. At first sight this is surprising, how- 
ever, we find that not the strength of the crossinterac- 
tion built in the crosspotential determines cluster mobil- 
ity. This is in accordance with recent findings in simi- 
lar systems for atomic deposition events 19 and for ion- 
beam intermixing 22 . The variation of parameter £ in Eq. 
(1) does not affect significantly the final results in accor- 
dance with earlier findings during ion- mixing of the Ti/Pt 
bilayer 22 . Parameter £ is proportional to the built in heat 
of mixing of the corresponding alloy phase 22 . Atomisti- 
cally, it determines the deepness of the potential energy 
well, hence the strength of the heteronuclear interaction. 
The possible reason of this insensitivity to the strength 
of the crossinteraction is that transient (ballistic) and 
athermal atomic transports are only weakly sensitive to 
pair-interactions e.g. between colliding ballistic particles 
e.g. in a collisional cascade. In other words, not chemi- 
cal interdiffusion (chemical forces) and thermochemistry 
(heat of mixing) determine the magnitude of bulk clus- 
ter mobility (and the implantation or penetration depth). 
CP occurs e.g. in Pt/Al either for weakly repulsive and 
for attractive crosspotentials. 



A. Mass-effect and cluster penetration 

It turned out during comparing the penetration depths 
obtained for various cluster/substrate pairs that atomic 
mass anisotropy, that is the atomic cluster to substrate 
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mass ratio (5) could seriously influence CP and could 
be responsible for the apparent asymmetry in CP with 
respect to the interchange of the cluster and substrate 
constituents. 

Since we find that atomic mass ratio could be an 
important ingredient of CP we briefly summarize our 
recent results in which we find the strong effect of 
mass anisotropy on interdiffusive behaviors of various 
film/substrate systems 22 . The enhancement of atomic 
penetration (intermixing) has also been reported recently 
of Pt in Ti substrate during ion-sputtering of Pt/Ti 
film/substrate system 20 . The asymmetry of interdiffu- 
sion in Pt/Ti with respect to the succession of the film 
and substrate has also been found 20 . In few other re- 
cent publication we also explained the enhanced atomic 
transport with mass-effect 22,3 . 

It might be the case that similar mass ratio driven ther- 
malization appears during CP. Indeed, if we set in artifi- 
cial mass isotropy in burrowing systems, we get a much 
weaker cluster penetration. Moreover, if we interchange 
artificially atomic masses between Al and Pt, we get no 
CP. Hence we find that the inversion of mass-anisotropy 
in the system hampers the enhancement of cluster mo- 
bility and CP. Moreover, if we interchange atomic masses 
in the Co/Au system, we does get, CP. In mass isotropic 
systems, in which no cluster penetration occurs, such as 
Ni/Cu or Co/Cu, ultrafast CP can also be induced by 
setting in artificial mass anisotropy. These simulation 
results strongly suggest that mass anisotropy could be a 
decisive parameter in cluster implantation. 

The role of mass effect in CP might be due to collisional 
cascade effects. When cluster impact occurs the initial- 
ization of local melting could not appear with atomic 



mass ratio S = T " ci < I, where m c i and m s 



arc the 



atomic masses in the cluster and in the substrate. This 
is simply due to backscattering effects when S < 1: light 
atoms stop at an interface composed of heavy atoms at 
low impact energies in a similar way as it has been found 
for intermixing during the ion-bombardment of various 
bilayers 22 . In the case of <5 ~ 1, however, the dissipa- 
tion of the impact energy into the lattice is the most 
effective 22 , and no thermalization or only a short lived 
collisional cascade takes place which does not allow CP. 
However, when 5 > 1, local thermalization of the sub- 
strate will appear because heavier cluster atoms pene- 
trate into the substrate top layer. 

The effect of S on CP can be demonstrated more quan- 
titatively if we vary 5 artificially e.g. in a nearly mass and 
atomic size isotropic system, such as Co/Cu and plot the 
penetration depth as a function of S. We use for simula- 
tions a standard Cleri-Rosato set of parameters 26 and a 
fitted potential for the Cu-Co interaction 36 . In principle 
we can get a universal plot, which holds more or less pre- 
cisely for any kind of cluster/substrate couples in which 
mass anisotropy governs CP. In Fig. 5 we show with 
closed symbols those values which are obtained for those 
couples, which composed of atoms with moderate lattice 
mismatch (e.g. in Cu/Ti, Pt/Al, Au/Al, Pt/Ag, Ti/Al 
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FIG. 5: The penetration depth (d pen in A) of various cluster 
impact events (2 eV/atom) as a function of mass ratio (6, 
5 > 1). Error bars denote standard deviations obtained dur- 
ing various events (statistics). Filled diamonds correspond 
to values obtained for the Co/Cu couple. The dotted line 
is for guiding eyes. Plus, circle, star, cross, triangle down 
and triangle up denote Cu/Ti, Pt/Al, Au/Al, Pt/Ag, Ti/Al 
and Au/Ag, respectively. Opened symbols: Results are shown 
with opened circle, diamond, triangle down, triangle up and 
square for Ni/Ti, Ni/Al, Co/Ti, Cu/Al and for Pt/Ti, respec- 
tively. A linear fit has also been plotted with a dashed line to 
guide eyes. Inset: The penetration depth (A) as a function 
of the lattice mismatch e (e > 0) for various cluster/substrate 
couples which are the followings: e = 0.024, Cu/Co, 0.058, 
Pt/Ti (iso), 0.04, Pt/Ag (iso), 0.10, Ni/Pd (iso), 0.11, Cu/Al 
(iso), 0.12, Cu/Ag (iso), 0.132, Cu/Ti, 0.14, Ni/Ag (iso), 0.15, 
Co/Ti, 0.15, Ni/Ti, 0.15, Co/Ti, 0.15, Ni/Ti, where (iso) de- 
notes artificial mass isotropic calculations. 



and Au/Ag). These systems give d pen values around the 
curve obtained for Co/Cu. 



B. The role of lattice missmatch 

However, in those cases, in which both cluster to sub- 
strate lattice mismatch (e = (a s — a c i)/a s > 0, where a s 
and a c i are lattice constants of the substrate and clus- 
ter constituents) and mass anisotropy are considerable 
(Co/Ti, Ni/Al, Pt/Ti, Cu/Al) a deviation from the curve 
of Co/Cu can be seen (opened symbols and dashed fit- 
ted curve). Hence we are faced with the splitting of the 
dpen vs. 5 plot into two regimes. These regimes can be 
understood as the separation of S and e dependent TCB 
processes. This is simply due to the fact that one hardly 
can find systems in which both S and e is robust. In 
strongly mass-anisotropic couples with 5 » 1 usually 
e is not that pronounced. When lattice mismatch is ro- 
bust in such a way that the cluster atom is the smaller 
(e < 0) mass anisotropy is not that robust as possible in 
other cases. Hence the frequency of the appearance of 
strongly mass anisotropic and lattice missmatched (<5, e) 
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couples among various possible cluster /substrate pairs is 
constrained by nature. 

In the nearly mass-isotropic Ni/Ti and Co/Ti only a 
smaller deviation can be seen, in Fig. 5 although still a 
considerable penetration occurs in these size-mismatched 
couples. In particular, the Co nanocluster burrows into 
the Ti phase temporarily completely and is ejected closer 
to the surface finally with d penet « 16 ± 1.5 A. However, 
no penetration takes place for the reversed cases (Ti/Ni, 
Ti/Co). 

We argue than in this paper that atomic size mismatch 
(ASM, also we call it atomic size anisotropy) could be 
together with 5 the basic parameters which govern cluster 
mobility and TCB in the bulk. 



C. Burowing clusters: atomic size mismatch and 
cluster bulk mobility 

Using MD simulations it has been shown recently 
that the ratio of the cluster lattice parameter to the 
substrate lattice constant (that is e) has significant 
effect on cluster diffusion on the surface (interfacial 
incommensurability) 2,25 . In the rest of the paper we 
argue that transient cluster burrowing and bulk mobil- 
ity could also depend on e. The possible operation of a 
e-dependent mechanism is supported by various simula- 
tions with mass isotropic and lattice mismatched systems 
(e.g. Co/Ti, Ni/Ti, etc.) which also show up transient 
cluster burrowing. 

Varying the elements in the cluster and substrate we 
can reach the conclusion that ASM could also be a key 
parameter in TCB. In particular, in the Ni/Cu system 
we get no TCB (e « 0.02), while if we replace Cu with 
Ti, TCB does occur in Ni/Ti (e w 0.15). The former cou- 
ple is a nearly mass and size isotropic system, while the 
latter one is size anisotropic being a considerable atomic 
volume difference between Ni and Ti. Also, no TCB oc- 
curs in the Ti/Ni system, when Ni is the substrate and 
Ti is the cluster. The larger Ti atoms can not pene- 
trate to the Ni phase. These situations support the ASM 
driven mechanism when S ss 1. Also, an artificially set 
mass isotropy e.g. in Pt/Al (e w 0.06) docs not sup- 
press completely CP, hence there must be another sys- 
tem parameter which drives CP even in mass-anisotropic 
couples. Co/Ti provides another mass isotropic and size 
anisotropic example in which TCB occurs. Finally, if we 
deposit Ti cluster on Al, no TCB occurs in accordance 
with the expectations (size isotropic couple). 

These results are summarized in the Inset Fig 5 in 
which the penetration depths of various cluster/substrate 
couples at 2 eV/atom impact energy have been plotted 
against the cluster to substrate lattice mismatch e. In the 
mass anisotropic cases in order to exclude mass effect on 
CP and to study purely the effect of lattice mismatch on 
CP an artificial mass isotropy has been imposed (iso). 
We assume that lattice mismatch and mass anisotropy 
have a nearly independent effect on CP, and their con- 



tributions to the total penetration depth d pen can be 
summed up from individual terms: d pen = d pen + d pen , 
where d pen and d pen are the mass anisotropy and atomic 
size mismatch induced contributions to transient cluster 
penetration. Hence during mass isotropic simulations 
dpen ~ d pen . We find a nearly linear dependence for 
various couples which supports the assumption that CP 
might depend on e. 

In few of the cases, however, deviation from the lin- 
ear dependence has also been found (Cu/Ag, Cu/Al). 
We find the strongest transient burrowing process in iso- 
Cu/Al while in iso-Cu/Ag the penetration depth drops 
below the fitted curve. In these cases it could be that not 
only anisotropy parameters govern CP, it could be that 
hidden (unknown) parameter(s) also contribute to CP. 
For vanishing lattice missmatch (e < 0.03) d p vanishes 
abruptly if 5 = 1 (mass isotropy). 



D. Nonburrowing clusters 

We summcrize hereby the general properties of non- 
burrowing cluster/substrate pairs. In general, in these 
couples cluster penetration and implantation can only be 
induced at least with few tens of eV/atom impact energy. 
For instance, in the Co/Cu system we find complete im- 
plantation at ~ 50 eV/atom energy, however, the cluster 
looses its integrity and even the substrate experiences se- 
rious damage. Hence, the technological applicability of 
this energy regime is questionable. 

As it has already been shown for many examples, the 
interchange of atoms in a burrowing cluster/substrate 
couple might lead to a nonburrowing couple. In these 
couples e < 0.0 (negative lattice mismatch) and/or 6 < 1 
(mass isotropy, 5 w 1 or inversed mass anisotropy, S < 1). 
The nonburrowing e < 0.0 couples are those A/B sys- 
tems in which A=A1, Ti, Ag, Au and B=Cu, Ni, Co, Pt. 
Hence e < 0.0 together with 5 < 1 have a strong suppres- 
sive effect on TCB independently from each other. The 
deposited cluster remains on the surface and d pen « 
when 5 < 1 and/or e < 0.0. In these cases penetration 
does not occur even temporarily. In many other mass- 
isotropic (<5 w 1) and size-anisotropic cases (e > 0.0) 
clusters first burrow deeply into the substrate (e.g. in 
Co/Ti) and than are ejected to the surface or to the top 
layers. 

The thermally induced burrowing of Co clusters into 
Au(lll) surface has been studied recently 24 . The ex- 
perimentalists claim that Co clusters possess a strong 
tendency to bury themselves into Au substrate 24 . Since 
this process can only be activated above room temper- 
ature annealing, we repeated our simulation for Co/Au 
with simulating the annealing process and indeed the par- 
tial burrowing of the Co cluster has been found at 450 
K, although the cluster does not keep its integrity com- 
pletely and sinks only few monolayer deep into the sub- 
strate. We conclude that in Co/Au the experimentally 
observed burrowing might not be an athermal process 
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(TCB) and could be similar to that of reported in the 
Co/Cu couple 13 . In the size-mismatched Co/Au system 
although e = 0.13, however, hence Co clusters are 

stopped at the Co/Au interface. The stopping power of 
S <C 1 is a stronger effect than the burrowing effect of 
lattice mismatch. Hence we conclude from this that the 
effect of ASM can only be realized when S > 1. This has 
been demonstrated with an artificially inverted S (the 
atomic masses of the cluster and substrate atoms have 
been interchanged) with which we can induce ultrafast 
CP in Co/Au. It is also clear that no penetration ap- 
pears in Au/Co, although 6 > 1. In this case this is due 
to the stopping power of e < 0. Hence S <C 1 and e < 
suppress TCB. 

V. CONCLUSION 

In conclusion, we explore the occurrence of tran- 
sient nanocluster penetration (burrowing) in various sub- 
strates using atomistic simulations at low energy clus- 
ter impacts. We point out that the cluster to substrate 
atomic mass and size anisotropy plays a significant role 
in the mass transport: clusters with heavier and smaller 
atoms burrow into a substrate composed of lighter and 
larger atoms. The transient burrowing process is largely 
insensitive to the strength of cross-interaction (that is 



proportional to the heat of mixing) and is driven by the 
atomic mass ratio and lattice mismatch. 

The deep penetration of metallic clusters into metal- 
lic substrates could allow the preparation of metallic in- 
clusions and burried nanostructures implanted into the 
substrate. Beyond the possible technological applica- 
tion of cluster burrowing, the explored new phenomenon 
could also be interesting in a theoretical point of view. 
In particular, the understanding of transient cluster mo- 
bility in the bulk could contribute to the advance of the 
emerging new field of anomalous diffusion (see refs. 19 ' 20 
and references therein). These peculiar results provide 
new evidences to the new concept that atomic trans- 
port could become anomalous in the nanoscale under not 
yet clearly established conditions 38 . Common feature of 
few of these processes the athermal characteristics (non- 
Arrhenius atomic transport) and that the deposit-surface 
interaction is largely independent of chemical forces 19,22 . 
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